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(54) Abstract Title 

Seismic detection apparatus with receivers defining a volume 

(57) Seismic detection apparatus comprising seismic detection means capable of detecting a plurality of 
seismic components over a defined tetrahedral volume (36) is provided. The seismic detection means 
comprises four three-component geophones (22, 24, 26, 32). Seismic data acquired by the geophones is 
processed to separate P-wave components from S-wave components. The geophones are spaced apart by 
distances smaller than the wavelength of the detected seismic components. The apparatus may be used on 
surface or in a marine environment or transition zone. A method of processing seismic data is also provided 
comprising acquiring seismic data relating to a wavefield over a selected volume of acquisition, and 
measuring the curl and divergence of the wavefield from the seismic data, to thereby identify seismic 
components within the seismic data. Additionally, an apparatus and method for hydrocarbon exploration is 
disclosed for using three or more seismic receivers placed in a plane and spaced closely to each other. 
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Seismic Detection Apparatus and Related Method 

This invention relates to seismic detection apparatus and to a 
5 method of analysing data acquired by such apparatus . 

BACKGROUND OF THE INVENTION 

In seismic detection, a seismic source signal propagates through 
10 different rock substrates or strata within the earth, so producing 
compression wave (P-waves) and shear wave (S-waves) energy which 
can be analysed to determine the direction and extent of geological 
features in strata. 

15 Generally a seismic source signal is produced, either on land, sea, 
or in a borehole, the source signal usually being produced as 
either acoustic or elastic energy. Some of the energy radiates 
downward through geological layers within the earth and is 
reflected in varying proportions at different layer boundaries and 

20 this reflected energy can be detected at the surface. Some of the 
energy within the seismic source signal propagates directly along 
the ground producing a wave signal known as ground-roll. To detect 
the acoustic reflections which are in the form of P-waves and S- 
waves, a linear or areal array of geophones is often used to 

2 5 reinforce the reflected energy and to attenuate the energy 

associated with ground-roll so that a good signal to noise ratio 
can be achieved for the reflected wave components of interest. 

It is desirable to be able to separate the P-wave and S-wave 

3 0 components of an elastic wave field as the separated components can 

be used more effectively to identify different characteristics than 
when combined. It is an aim of the present invention to provide 
improved separation of seismic signals into P-V7ave and S-V7aves . 



Methods have been described which use three or more sensors in 
close proximity (e.g. less than 1 meter) for determining certain 
ground characteristics very close to the sensors. For example, in 
5 published European Patent Application 0 455 091 A2 (Application No. 
91106459.0)/ a method is described which uses relatively closely 
spaced '^oscillation sensors" to determine the "laminar structure 
and other characteristics of [the] ground." However, this method 
only provides info2rmation about characteristics of the ground 

10 within a distance from the sensor array of the order of the maximum 
wavelength of the surface waves (hereinafter referred to as the 
''near-site grouind structure"). Moreover, the method relies on 
sensing "microseisms" which is believed to refer to naturally 
occurring seismic activity. Thus, the published application does 

15 not disclose a method using seismic energy from a man-made seismic 
source, nor does it teach a method of sensing the components used 
to process reflected waves . Rather the method is only concerned 
with groimd-roll type waves and does not teach or suggest a method 
suitable for hydrocarbon exploration using seismic energy reflected 

20 from man-made seismic sources. 

SUMMARY OF THE INVENTION 

2 5 According to one aspect of the present invention, there is provided 

seismic detection apparatus comprising seismic detection means 
capable of detecting a plurality of seismic components over a 
selected volume. 

3 0 Preferably the seismic detection means comprises a plurality of 

receivers spaced apart to enclose the selected volume. The 
receivers may include hydrophones, geophones, accelerometers or a 
combination of these, and preferably are adapted to detect seismic 



- 3 - 



signals along three mutually orthogonal axes. Thus the receivers 
may be provided by three- component geophones , 

Preferably the selected volume is formed by at least three 
5 receivers spaced within a plane and at least one receiver placed 
external to the plane, thereby to enclose the selected volume over 
which measurement of seismic signals occurs . 

In accordance with the invention, the receivers can be spaced to 
10 form a tetrahedron with three receivers spaced to define a triangle 
within the plane and a fourth receiver placed external the plane 
but in line with the centroid of the triangle. Other 
configurations of spaced apart receivers may be used depending on 
the characteristics of the wavefield being measured, for example 
15 receivers may be spaced to form an octahedral volume, a cubic 

vol-ume, or a spherical voliime. Certain shaped volumes may be more 
appropriate for different uses, for example receivers spaced to 
form an octahedral volume may be more appropriate for use with 
permanent downhole sensors . 

20 

The seismic detection apparatus in accordance with the present 
invention is suitable for use on the ground surface or ocean 
bottom. When used on the ground surface, the enclosed volume can 
be achieved by placing one receiver further below the surface of 

2 5 the other receivers, for example by burying it beneath the surface. 
With marine and ocean bottom surveys, the volume may be defined by 
tethering the receivers at different water depths, for example by 
placing three receivers spaced apart at the same depth, and a 
fourth receiver at a different depth. The seismic detection 

30 apparatus according to the present invention is also of use in 

seismic detection in transition zones such as swamps and marshes, 
where again tethering of the receivers at different depths will 
define the enclosed vol-ume. 
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The apparatus may also have all the receivers incorporated into one 
body for attaching to a wireline and using downhole, or for using 
permanently downhole. Such a permanent sensor is of particular 
5 advantage where monitoring of a reseirvoir or other geological 
feature over time is required. This may be appropriate when 
monitoring fluid flow occurring within a reservoir as permanent 
sensors can be used to detect natural seismicity resulting from 
fluid opening and closing micro-fissures within the substrates and 

10 producing micro-earthquakes. Whilst these micro-earthquakes 

produce very weak broadband signals, by use of permanent sensors, 
this activity can be detected and production adapted by remedial 
action such as acid to remove carbonate deposits, or fracturing 
rock with high pressure steam or water, so as to improve production 

15 efficiency. 

Preferably the apparatus further comprises processing means which 
analyses detected seismic components from individual receivers to 
separate P-wave components from S-wave components. The processing 

2 0 means may be provided at substantially the same location as the 

receivers to allow for on-site processing, or remote from the 
receivers to allow off -site processing. 

In surveys of the vertical component of seismic data, the receivers 
25 may be adapted to detect components along one axis and orientation 
of the receivers selected so as to permit processing of one- 
component data to identify P-wave components. Therefore the 
receivers may be provided by one-component geophones. 

3 0 Preferably the spacing of the receivers is selected to be smaller 

than the wavelength of the detected seismic components. A 
preferred spacing is less than 1 m, more preferably between 0.05m 
and 0-50m, although larger spacings of up to a third of the 



shortest detected wavelength are possible, and thus the range may 
be between 1.5 and ISm. 

In use, the apparatus may be operated to detect seismic signals 
5 over a time period of the order of five seconds where P-wave 
components are of principal interest, with any marginal S-wave 
component detected at that time being used to produce improvements 
in the signal to noise ratio. 

10 Where a detection time of the order 10 to 15 seconds is used, both 
p_;^ave and S-wave components will be detected within the seismic 
signals and may be separated for multi-component inversion yielding 
images of reservoirs and other geological formations, and full 
reservoir characterisation . 

15 

In accordance with another aspect of the invention, there is 
provided a method of processing seismic data comprising acquiring 
seismic data relating to a wavefield over a selected volume of 
acquisition, and measuring the curl and divergence of the wavefield 
2 0 from the seismic data, to thereby identify seismic components 
within the seismic data. 

As p-waves are curl free and S-waves divergence free, this allows 
the separate identification of P-wave and S-wave components. Up- 
2 5 going and down-going wavefield components in boreholes may also be 
identified from the seismic data. Attenuation of unwanted seismic 
components from data acquired downhole is also enabled, and noise 
sources such as ground roll, airwave and inter-sensor perturbations 
can be removed from the data by signal processing. 

30 

The method preferably further comprises averaging the curl and 
divergence over the selected voliame of acquisition, and where 
appropriate acquiring seism.ic data over a selected period of time. 
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Additionally the method may further comprise defining the selected 
voltmie of acqiaisition to be small compared to the wavelength of 
seismic waves to be detected. 

5 

As described above, in a seismic acquisition at a free surface, the 
constraints imposed by the free surface enable the reduction of the 
acquisition volume to a acquisition surface where all receivers 
are located in a plane. 

10 

In another aspect of the invention, the representation theorem 
based filters used for removing multiples from Ocean Bottom Cable 
data as described in the Published UK Patent Application GB 2 33 3 
364 A (Published 21 July 1999) is modified to provide a method of 
15 separating the seismic data acquired with the planar receiver 
configuration above into up- and downgoing wavefield. 

In a variant that is particularly suitable for hydrocarbon 
exploration seismic acquisitions at a " free surface or fluid/solid 

20 boundary, the volume formed by the receivers can be reduced to a 
plane. Hence, three 3 -component receivers can be used in 
combination with the free surface (or fluid/solid) condition 
constraints to advantageously determine approximations to the 
spatial and temporal derivatives local to the receivers 

25 (hereinafter the "full wavefield"). The receivers should be 

adapted to measure the seismic energy from a seismic source that 
has been reflected from one or more subterranean earth structures 
outside the near-site ground structure. 

30 According to a further aspect of the invention, each of the 

receivers are one- component geophones and a fourth seismic receiver 
is located less than about 1 meter from the other three seismic 
receivers . 
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A method for hydrocarbon exploration is also provided that includes 
imparting seismic energy into the earth that travels through the 
earth and reflects off of one or more subterranean earth structures 
5 outside the near-site ground structure. At least three seismic 
receivers are positioned substantially in a plane at or near the 
earth surface, and the seismic receivers are spaced less than about 
1 meter from each other to form a receiver cluster. The received 
seismic energy is stored and analysed to determine characteristics 
10 of one or more hydrocarbon reservoirs. 

Preferably, the curl and divergence of the wave field is determined 
from the seismic data, thereby identifying seismic components 
within the seismic data. 

15 

Preferably, the P-wave and S-wave components are separately 
identified, the up-going and down-going wavefield components are 
identified from the seismic data, and unwanted seismic components 
from seismic data is attenuated. 

20 

Preferably, the curl and divergence is averaged over the selected 
plane or volume of acquisition, and the selected plane or volume of 
acquisition is defined to be small compared to the wavelength of 
seismic waves to be detected. 

25 

These and other features of the invention, preferred embodiments 
and variants thereof, possible applications and advantages will 
become appreciated and understood by those skilled in the art from 
the detailed description and drawings following hereinafter. 

30 

BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1 shows an example of the distribution of seismic data over 
a specified period of time; 

FIG. 2 shows an example of a seismic data acquisition and system 
5 according to the present invention; 

FIG. 3 shows one embodiment of seismic detection apparatus 
according to the present invention; 

10 FIG. 4 shows a second embodiment of seismic detection apparatus 
according to the present invention; 

FIG. 5 shows a tetrahedron used in explanation of a processing 
method according to the present invention; 

15 

FIG. 6 shows a planar configuration of receivers according to 
the present invention; 

FIG. 7 shows a top view of a configuration of four IC geophones 

2 0 in a plane that can be used to measure the divergence on the 

Earth's free surface according to the present invention;, and 

FIG. 8 summarises a processing method based on a planar receiver 
arrangement according to the present invention. 

25 

MODECS) FOR CARRYING OUT THE INVENTION 

In Figure 1, a typical hodogram of a three component marine shot 
gather survey is shown, where the x axis is the offset in meters 

3 0 and the y axis is time in seconds from the generation of a source 

signal at t=0. A hodogram displays the different seismic 
components detected over a specified time-window. The first 
arrivals around t=ls predominantly comprise P-waves 10, 10'. Near- 
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surface scattering and channelling in slower near-surface layers 
are revealed in the P-wave arrivals contained in the hodogram 
through some degree of elliptical motion. Following the P-wave 
arrivals is a relatively quiet section 12, 12' (where deeper 
5 reflections can be observed) until the S-wave arrivals 14, 14' are 
observed near t=3s. The S-waves also scatter and become channelled 
in the near-surface which causes ellipticity in the hodogram. 

Following the S-wave arrivals are elliptically polarised Rayleigh 
10 waves 16, 16' (ground-roll). The innermost ''cone" 20 of the 

section consists of strongly scattered energy that typically masks 
deeper reflections . 

Noise in seismic data to a large degree comprises source-generated 
15 ground-roll and S-waves. Separating P-wave data from S-wave data 
is useful to reduce noise attenuation. Also P-wave data may be of 
primary interest with the complete removal of S-waves being 
preferred, for example in shallow seismic imaging for engineering 
or environmental purposes . 

20 

Measuring the P-wave and S-wave components of an elastic wavefield 
separately has several advantages. First, separation of P-waves 
and S-waves enables better and more accurate analysis of each 
separated component thus improving moveout analysis and Amplitude 
25 Versus Offset (AVO) analysis. Secondly heterogeneity in the near- 
surface (as well as deeper down in the crust) and rough surface 
topography, scatter S-waves to a greater extent than P-waves, with 
scattered S-waves significantly contributing to noise in seismic 
reflection data. 

30 

Figure 2 is a schematic illustration of a seismic data acquisition 
and processing system. Seismic sources 110, 112, and 114 are 
depicted which impart seismic energy into the earth at its surface 
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126 at predetermined times upon command of a controller (not 
shown) . The seismic sources may comprise impulsive sources (such as 
explosives, or in the marine environment, airguns) or vibratory 
sources. While three sources are depicted for simplicity, other 
numbers of sources, including a single source, could be used. The 
seismic data acquisition could be carried out on land, marine, or 
transition-zone environments. In the case of marine acquisition, 
12 6 would depict the ocean bottom, and the sources 110, 112 and 114 
would typically be submerged at a predetermined depth near the 
water surface (not shown) . In this case receivers 116, 118, and 12 0 
can be in an ocean bottom cable. In the case of land acquisition, 
surface 12 6 is the ground, and truck mounted vibrators can be used 
as sources 110, 112 and 114. The seismic energy imparted by the 
sources travels through the earth in the form of both p and s 
waves; this is schematically depicted in Figure 2 as, for example, 
arrow 134. The energy reflects off of certain subterranean 
surfaces, here depicted as surface 122 and surface 124, and 
eventually reach and are detected by receivers 116, 118, and 12 0 as 
depicted, for example, by arrow 138, According to the preferred 
embodiment, receivers 116, 118 and 120 will be geophones, but other 
suitable receivers such as optical sensors could be used. Note 
that receivers 116, 118 and 12 0 in Figure 2 correspond to the 
receivers and geophones depicted in Figures 3, 4, 6, and 7. 
Although only three receivers are shown in Figure 2, there could be 
more receivers used in a particular embodiment- For example, in 
the embodiments shown in Figures 3, 4 and 7, at least four 
receivers are used. 

Referring again to Figure 2, the apparatus shown is being used for 
seismic exploration of hydrocarbon reservoirs. Cap rock 12 8 is 
relatively impermeable to the hydrocarbon substance. The 
hydrocarbon substance is thereby trapped in hydrocarbon reservoir 
130. Hydrocarbon reservoir 130 is typically made up of a reservoir 
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rock (for example sandstone in a sand channel) and the interface 
between it and the cap rock 12 8 reflects the seismic energy back to 
the receivers as depicted by arrow 138. Also shown in Figure 2 is 
ground roll (e.g. Rayleigh waves) which is schematically depicted 
5 as arrow 13 6. Ground roll travels in the near-surface region and 
does not penetrate to depths necessary for locating and 
characterising hydrocarbon reservoirs such as reservoir 130. 

Receivers 116, 118, and 12 0 convert the seismic energy into 
10 electrical signals and transmit these signals to a central 

recording unit 140. In the case of land acqpaisition this central 
recording unit 140 is typically located at the local field site, 
and in marine acquisition it is typically located on the vessel 
used to tow the sources and receivers (not shown) . Once the central 
15 processing unit 140 performs correlation and other desired 

processing, it typically stores the data in the form of time-domain 
traces on a magnetic tape. The data, in the form of magnetic tape 
is later sent for processing and analysis to a seismic data 
processing centre, typically located in some other geographical 
20 location. The data transfer from the central recording unit 170 in 
Figure 2 is depicted as arrow 176 to a data processor 180. 

Figure 3 shows one embodiment of a seismic detection apparatus 
according to the present invention, with four three- component 

2 5 detectors, such as geophones or hydrophones, being placed to form 
an enclosed tetrahedral volume. Geophones are velocity motion 
sensors, with hydrophones being pressure sensors. Geophones 22, 
24, 2 6 are arranged in a horizontal plane 3 0 with each geophone 
being spaced the same distance from the other two geophones, so as 

30 to form an equilateral triangle. The fourth geophone 32 is placed 
vertically below the centroid 34 of the equilateral triangle formed 
by geophones 22, 24, 26, Placing of geophone 32 external to the 
horizontal plane 3 0 is achieved by positioning the horizontal plane 
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30 on or close to the Earth's surface with the fourth geophone 32 
dug to a position well below surface. Within a marine, ocean 
bottom or swamp environment, the arrangement of receivers is 
achieved by tethering three hydrophones and/or geophones at one 
5 depth and one hydrophone and/or geophone at a different depth. 

• An ideal tetrahedral vol-ume is not essential, and thus for example 
two geophones placed external to horizontal plane 3 0 can be used. 

10 The geophones 22, 24, 2 6 are positioned one meter apart, although 
spacing will vary depending on the surrounding rock substrate and 
how fast waves travel within the substrate. The spacing of the 
geophones to form an enclosed volume 36 is chosen so that the 
wavefield is detected on a scale which is smaller than the 

15 wavelengths of interest. For wave freq[uencies of between 50 to 10 0 
Hz with velocity ranging from 300-2000 m/s, the wavelength tends to 
be in the range 5 to 50 m. To record the wavefield over such a 
range, the spacing of the geophones is generally a third of the 
wavelength and generally in the range 1.0m to 15m, although may be 

20 as small as 0.05 m. The tetrahedron arrangement shown in Figure 3 
uses the minimum n\imber of geophones required to define an enclosed 
volume, although other configurations with more geophones are 
possible, such as those forming an octahedron,- cube or sphere. 

25 The use of four three-component geophones distributed as shown in 
Figure 3 allows the separation of the wavefield into P-waves and S- 
waves . After the seismic data has been recorded either where 
acquisition of data is occurring or remote from the geophones, 
differences in geophone coupling can be removed by normalising 

30 separately the responses of each component at the corners of the 

tetrahedron. If required, more accurate compensation for coupling 
differences can be obtained by using information from adjacent 
tetrahedral arrays of geophones, while preserving the offset 
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varying character of the seismic data. The geophone coupling may 
also be frequency dependent varying from geophone to geophone. 

The three- component geophones 22, 24, 26, 32 detect components of 
5 the seismic signal at three orthogonal axes. By measuring the curl 
of a displacement field u (x, t) within the enclosed volume 3 6 and 
the divergence of the displacement field within the volume 36, it 
is possible to separate out the P-wave and S-wave components as P- 
waves have zero curl and S-waves have zero divergence. 

10 

With the array of geophones shown in Figure 4, the four geophones 
40, 42, 44, 46 are one-component geophones. The array encloses a 
tetrahedral volume 5 0 once more, but the axis along which each one- 
component geophone measures is calculated and precisely implemented 

15 when positioning the apparatus to ensure that the divergence can be 
deduced from the single components measured. Measuring with single 
component geophones provides less information about the seismic 
signals. However at present 90% of seismic surveys only monitor 
the vertical component, and the apparatus of Figure 4 when used 

2 0 with a listening time of five seconds allows a substantial 

improvement in the signal to noise ratio if looking principally for 
P-wave components. 

Once data are acquired using the geophone array enclosing a volume, 
2 5 separation of the wave- field into ?-waves and S-waves or 

improvement of the signal to noise ratio for P-wave components, can 
be undertaken by use of the method described below. 

The elastic wave equation for a displacement field M{x,t} in an 
30 isotropic medium with Lame parameters A and and density p can be 
written as: 



[1] 



p U = f + (X -f 2fX) V(V • u) + ^iV X (V X u) 
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Here f denotes a distribution of body forces. Lame's thereon (Aki 
and Richards, 1980) states that there exist potentials ^ and Y of 
u with the following properties: 



[2] u = V4> + V X T , 

[3] V • Y = 0, 

[4] ^ = — -h a^V^O , 

P 

[5] ^ — -h pV^Y, 

p 



where a and p are the P- and S- velocities. An elastic wavefield u 
can thus be decomposed into its P- and S-wave components, V<E> and V 
X respectively. Equations [2] and [3] yield 

[6] V • u V^O , 

[7] Vxu = VxVxY 

(R&de, L., and Westergren, N. , 1988, BETA p Mathematics Handbook: 
Student litteratur, Lund.). By measuring the curl and the 
divergence of an elastic wavefield we can thus measure the P- and 
S-wave components separately. 

The technique for P- and S-wave separation is based on recording 
the elastic wavefield around a small volume as compared to the 
wavelength. The volume may be of arbitrary shape, but the 
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equations for wavefield separation are simplified if the volume has 
a known deterministic shape such as that of a tetrahedron, cube, 
octahedron or a sphere. A co-ordinate system that is appropriate 
for the geometry employed simplifies the equations further. The 
5 discussion of the method below uses a tetrahedron and the formulae 
are derived in Cartesian co-ordinates . 

Figure 5 shows a tetrahedron with corner points Pi, P3 / and P4 . 

We denote the Cartesian co-ordinates of point Pi by (xi, y^, zi) . 

10 Let PjP^ denote the vector that starts at Pi and ends at Pj . The 

volume of the tetrahedron in Figure 5, V^ecra/ can be calculated from 
the absolute value of the determinant of the matrix M with the 

vectors P2P1 / P3P1 / and P4P1 as rows: 

15 [8] Vtetra = " pet M| . 

5 

Gauss' theorem states 

[9] JJj;(v.u)dv = J£u.dS 

20 or alternatively. 



[10] JJI;^^ X u)dV = j£u X dS 



If the dimensions of the tetrahedron configuration in Figure 5 are 
25 small compared to the wavelength of u, the mean values of the curl 
and divergence in equations [9] and [10] can be calculated from: 



[11] V . u = -r^|JJ,(V - u)dv 



3 0 and 
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10 



[12] V X u = — — lllj^ ^ u)dV 



^tetra 



Moreover, the surface integrals on the right hand side of equations 
5 [9] and [10] reduce to four surface integrals respectively, over 
plane triangular surfaces . 

Since the normal to the surface of the tetrahedron does not vary- 
over each of the triangular sub-surfaces, equation [9] becomes 



[13] //j;(v .u)dv = j£_^^udsj.,,3 



u dS • n243 

^243 



15 + JJ u dS • ni42 



■'Si 4 2 



+ J[ u ■ ni34 , 



•'Si 3 4 

where nijjt is the outward unit normal to the triangular subsurface 
2 0 S±j}c with corner points P^, Vj and Pjt (see Figure 5) . 

Where the side of the tetrahedron is small compared to the 
wavelength of the elastic wavefield, the value of u at any point on 
one of the triangular surfaces of the tetrahedron can be calculated 
25 by linear interpolation of the wavefield between the three corner 
points of the triangle. The surface integral over Sijk (using 
linear interpolation between P^, P-,- and P;c) then reduces to 
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[14] ff udS = ^ (x^ + Xj +Xk, Yi + Yj + Yk' Zi + + Zj,), 



where A is the area of Sijk- Since Aiii-jk = P-iPi ^ P^Pi /2 we have 



5 [15] JJyV • U) dV = - [ 



(xi + X2 + X3, yi + y2 + y3. -zi + 22 + ^3) (^2^1 ^ P3P1 ^ 



+ {X2 + Xi + X3, y2 + y4 + y-i. zz + ^4 + zs) ( P4P2 ^ P3P2 ) 

10 



+ (xi + X4 + jf2, yx + yi + y2, zi -s- Z4 Z2) ( P4P1 x P2P1 ) 



+ {xi + X2 + Xi, yi + y3 + yi, z-, + Z3 + Z4) (P3P1 x P4P1 ) ] 
15 Equivalently, the curl of u in V can be calculated as 

[16] JJJJ^(Vxu)dV = ^[ 



20 
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Finally, the average values of the divergence and curl of u within 
V are obtained from equations [11] and [12] . 

In this way, the divergence and curl of a wavefield can be 
5 calculated from a set of closely spaced measurements distributed in 
all spatial directions. A minimum of four measurement points is 
needed. The volume that is bounded by planes connecting these 
points is generally a tetrahedron and where the dimensions of the 
tetrahedron are small with respect to the wavelength of the 
10 wavefield, an assumption that the wavefield varies linearly along 
each side of the enclosed volume is possible. The average 
divergence and average curl can then be efficiently calculated 
using Gauss' theorem. 

15 For instance the average divergence inside the tetrahedron is 
calculated from the volume integral of the divergence of the 
wavefield over the tetrahedron, divided by its volume. The volume 
integral can in turn be converted to four surface integrals over 
each side of the tetrahedron by using Gauss' theorem. If the 

2 0 wavefield varies linearly over each side, these surface integrals 

are simply the volume of a prism where the base of the prism is the 
surface area and the heights of the three corners are the magnitude 
of the wavefield in the respective corners projected onto the 
normal of the surface . 

25 

The average curl of the wavefield is evaluated similarly. The same 
procedure can be used to calculate average curl and divergence over 
more general volumes. 

3 0 Uncertainties in the actual geometry can be eliminated by allowing 

the exact locations of Pi, P2/ P3/ and P4 to vary subject to either 
minimising equation [16] over a window(s) with an identified P- 
event and/or minimising equation [15] over S-events. This 
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technique may also serve as an alternative means for adjusting for 
the differences in geophone coupling over the tetrahedron. 

For anisotropic media the decomposition of u into P- and S-wave 
5 components as described by equations [2], [3]^ [4] and [5] is 

strictly not valid since the (quasi) P- and S-waves in anisotropic 
media both have compressional and shear components. However, this 
is of secondary importance and the theory above will be 
approximately valid for most practical cases. 

10 

In heterogeneous media, the heterogeneities act as secondary 
sources and scatter the wavefield. The separation of a wavefield 
into P-waves and S-waves is not strictly valid in the vicinity of a 
source. However the P/S separation technique described herein is 
15 generally applicable for common heterogeneous and anisotropic Earth 
materials, as well as for isotropic materials as discussed above. 

With the apparatus and associated method, P-wave and S-wave 
components in seismic data can be separated. 

20 

In addition attenuation of near-surface scattering is possible as 
much noise is of S-wave character and where the primary interest is 
to record P-wave data, S-waves can be removed by employing the 
configuration of four one-component geophones as shown in Figure 4. 

25 

The present invention is particularly well-suited for seismic 
techniques applied to hydrocarbon exploration. However, the 
present invention also finds application in environmental, 
engineering and groundwater exploration and imaging. For such 
30 applications, it is typically of interest to image or detect very 
shallow structures, as shallow as 1-lOOm. Near-surface generated 
noise is particularly severe in these applications and a 
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configuration of receivers such as in Figure 4 is particularly 
advantageous for shallow seismics. 

The invention also allows for the attenuation of ground-roll. 
Since ground-roll can be elliptically polarised, it will appear 
both in the separated P- and S-wave sections, but with lower 
amplitude in the separated P-wave sections. 

P- and S-wave separation is of particular use for marine seismic 
data recorded at the seafloor, since the near-surface at the 
seafloor is generally is of a more homogeneous character than at 
land. The technique is of use in removing Scholte waves that 
typically contaminate ocean data. Moreover, measurements of the 
wavefield over a volume also allows for the separation of the 
wavefield into up- and down-going waves . 

Figure 6 shows a receiver arrangement particularly useful for land 
acquisition. It comprises three 3C (component) geophones 52, 54, 56 
arranged in a plane on the surface of the survey terrain. Each 
component of each geophone is measured independently. Thus, the 
three geophones are connected to the recording or processing unit 
58 via a total of nine channels 60. Cable 62 provides a signal 
transmission between the groups of geophones. 

The preferred spacing of the geophones in this embodiment is chosen 
in the same manner as described above, i.e. so that the wavefield 
is detected on a scale which is smaller than the wavelengths of 
interest. For wave frequencies of between 50 to 100 Hz with 
velocity ranging from 3 00-2000 m/s, the wavelength tends to be in 
the range 5 to 50 m. To record the wavefield over such a range, 
the spacing of the geophones is generally a third of the wavelength 
and generally in the range 1.0m to 15m, although may be as small as 
0 . 05 m. 



) 
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In this embodiment, the calculations according to the above 
eqriations are modified by introducing the free surface boundary 
condition. This condition constrains the possible solutions for the 
5 wavefield such that it is not necessary to employ a (fourth) out- 
of -plane receiver to sample additional components of the wavefield. 

Furthermore the free surface condition is introduced into a 
formalism known as the elastodynamic representation theorem. This 
10 theorem is used for example in the Published UK Patent Application 
GB 2 333 3 64 A (Published 21 July 1999) to derive filters for 
separating up- and downgoing wavefield in the case of ocean bottom 
acquisition. 

15 Using the representation theorem and the boundary conditions of a 
solid-air interface, the following equations for up/down going 
separation of particle velocities v recorded by the planar receiver 
arrangement can be obtained: 

2 0 [17] v:^(x) = - ( v^(x) - f;;;(x) * V3(x) ) , 

[18] v5(x) = - { V3(X) - RT^X) * Vy(x) ) . 

The subscripts n and m denote indices for the horizontal coordinate 
25 axes x and y (for the 2-D case they reduce to x) . The filters in 
equations [17] and [18] are 

[19] f;;:^ = -2B^{i + 2k^^dl) go^^^f + 2kp2a2g(kp) and 



3o' [20] f;;:^ = 2]cp2a^(_2a2g(]c„) + oc| + 2di)go^p) 
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In particular for the two dimensional (2D) case ( in the frequency- 
waveniimber domain), using expressions for Green's functions in a 
homogeneous medium , equations [19] and [20] become: 



[21] 
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where — -yk^^^ — k^^ is the vertical wavemamber for P-waves and 

k^^' = ^kp — k^ is the vertical wavenumber for S-waves. 



Filters for combinations of partial derivatives of velocities (e.g. 
divergence and curl) can be obtained using equations [17] and [18]. 
Note that for an isotropic elastic mediaxm vertical derivatives of 

15 particle velocities recorded at the free surface can be converted 
to horizontal derivatives (i.e. the * f ree surface condition' , see 
J. 0. A. Robertsson, ''A numerical free-surface condition for 
elastic/viscoelastic finite-difference modeling in the presence of 
topography". Geophysics, 1996, Vol. 61, pp. 1921-1934, incorporated 

2 0 herein by reference) . A similar transformation of vertical to 
horizontal derivatives can be performed at the ocean sea bottom 
using the fluid/solid boundary condition. A related approach has 
been disclosed in the Published UK Patent Application GB 2 333 364 
A (Published 21 July 1999) , incorporated herein by reference. 

25 

According to another embodiment, shown in Figure 7, one-component 
sensors can also be used in a closely spaced surface-plane 
arrangement. Figure 7 is a top view of a configuration of four IC 
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geophones, 80, 82, 84 and 86 in a plane that can be used to measure 
the divergence at the Earth's free surface. As mentioned above 
vertical derivatives can be converted to horizontal derivatives 
using the free surface condition (or in a marine environment, 
5 fluid/ solid boundary condition) . Therefore, by combining single 
components of horizontal particle velocity measured by different 
geophones the divergence can be calculated. Similar configurations 
of IC geophones can be designed to measure curl components and/or 
to include up/down separation of divergence and curl. 

10 

Additionally, the present invention can be implemented as a cluster 
of sources, rather than as clusters of receivers, according to 
theorem of reciprocity. In this embodiment the spacing and 
relative placement of the receivers in Figures 3-6 would instead 
15 apply to the spacing and placement of the sources. The analytical 
method described above could then be used to achieve the objects of 
the invention. 

While preferred embodiments of the invention have been described, 
2 0 the descriptions are merely illustrative and are not intended to 

limit the present invention. For example, although the embodiments 
have been described in the context of using geophones as the 
receivers of seismic energy, other types of receivers are 
contemplated and could be used in the present invention. For 
25 example in Figure 6, receivers 52, 54 and 56 could of the type 
based on optical and/or laser beams. 



30 
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CLAIMS 



5 1, A Seismic detection apparatus comprising: 

a first, second and third receiver, each of the receivers 
spaced within a plane; and 

a fourth receiver located external to the plane thereby 
enclosing a selected volume bounded by the first, second, 
10 third and fourth receivers over which measurement of seismic 

signals occurs . 

2. A seismic detection apparatus according to claim 1 
wherein the receivers comprise geophones . 

15 

3. A seismic detection apparatus according to claim 1 
wherein the receivers comprise hydrophones , 

4. A seismic detection apparatus according to claim 2, 
20 wherein the receivers comprise accelerometers . 

5. A seismic detection apparatus according to claim 2, 
wherein the seismic detector is adapted to detect seismic signals 
along three mutually orthogonal axes . 

25 

6. A seismic detection apparatus according to claim 1, 
wherein the receivers are spaced to form a tetrahedron, with three 
receivers spaced to define a triangle within the plane and a fourth 
receiver placed external the plane but in line with the centroid of 

30 the triangle. 



7. A seismic detection apparatus according to claim 2, 
further comprising a processor adapted and configured to analyse 
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detected seismic components from individual receivers to separate 
P-wave components from S-wave components . 



5 wherein the processor is provided at substantially the same 

location as the receivers so as to allow for on-site processing. 



wherein the processor is provided in a location remote from the 
10 receiver so as to allow off-site processing of detected seismic 
components . 

10. A seismic detection apparatus according to claim 1, 
wherein the receivers are adapted to detect conroonents along one 

15 axis and orientation of the receivers is selected so as to permit 
processing of one-component data to identify P-wave components 
within the plurality of seismic components. 

11. A seismic detection apparatus according to claim 1, 

2 0 wherein the spacing of the receivers is selected to be smaller than 
the wavelength of the detected seismic components. 



8. 



A seismic detection apparatus according to claim 1, 



9. 



A seismic detection apparatus according to claim 7, 



12. A seismic detection apparatus according to claim 1, when 
attached to a wireline for placing temporarily downhole. 



25 



13. An apparatus for hydrocarbon exploration comprising: 



a seismic source adapted to impart seismic energy into 



the earth at a predetermined time; 



at least three seismic receivers each adapted to 



30 



measuring the seismic energy from the seismic source that has 
been reflected from one or more subterranean earth structures 
outside the near-site structure, the seismic receivers located 
substantially in a plane at or near the earth surface, v/herein 



said at least three seismic receivers are spaced less than 
about 1 meter from- each other to form a receiver cluster , 

14. An apparatus according to claim 13 further comprising a 
fourth seismic receiver located less than about 1 meter from the at 
least three seismic receivers and wherein each of the receivers are 
one -component geophones . 

15. An apparatus according to claim 14 wherein the each of 
the one-component geophones are positioned in a diamond- shaped 
pattern and each of the sensed components are in different 
directions within the plane. 

16. An apparatus according to claim 13 wherein each of the 
receivers are geophones capable of measuring at least three 
components . 

17. An apparatus according to claim 16 wherein the earth 
surface is the surface of land. 

18. An apparatus according to claim 16 wherein the earth 
surface is the sea bottom. 

19. An apparatus according to claim 13 wherein a hydrocarbon 
reservoir is located near the one or more subterranean earth 
structures outside the near-site structure. 

20. An apparatus according to claim 16^ wherein the at least 
three components of each of the geophones are transmitted to a 
processing unit to calculate an approximation of the full wavefield 
at the position of the cluster making use of a free surface 
bo-undary condition. 
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21, A method for hydrocarbon exploration comprising: 
imparting seismic energy into the earth at a predetermined 

time such that the seismic energy travels through the earth and 
reflects off of one or more subterranean earth structures outside 
5 the near-site structure; 

positioning at least three seismic receivers 
substantially in a plane at or near the earth surface, wherein 
the seismic receivers are spaced less than about 1 meter from 
each other to form a receiver cluster; 
10 receiving seismic energy imparted by the seismic source 

that has been reflected off of the one or more subterranean 
earth s true tur es ; 

storing seismic data that represents the received seismic 
energy; and 

15 analysing the data such that characteristics of one or 

more hydrocarbon reservoirs can be determined. 

22, A method according to claim 21 wherein the step of 
analysing comprises the step of measuring the curl and divergence 

2 0 of the wavefield from the seismic data, thereby identifying seismic 
components within the seismic data. 

23, A method according to claim 22, V7herein P-wave and S-wave 
components are separately identified. 



25 



24, A method according to claim 23, wherein up-going and 
down-going wavefield components are identified from the seismic 
data. 
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25. A method according to claim 24, wherein the step of 
analysing further comprises attenuating unwanted seismic components 
from seismic data. 
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26- A method according to claim 25, wherein the step of 
analysing further comprises the step of averaging the curl and 
divergence over the selected plane or volume of acquisition. 

5 27. A method according to claim 26, wherein the step of 

analysing further comprises the step of defining the selected plane 
or volume of acquisition to be small compared to the wavelength of 
seismic waves to be detected. 
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